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Abstract. The investigation of lighit nuclei with ab-initio methods provides an optimal setting to probe our 
knowledge on nuclear forces, because the few-nucleon problem can be solved accurately. Nucleons interact not 
only in pairs but also via many-body forces. Theoretical efforts need to be taken towards the identification of 
nuclear observables sensitive to the less known many-nucleon forces. Electromagnetic reactions can potentially 
provide useful information on this. We present results on photo-absorption and electron scattering off light nu- 
clei, emphasizing the role of three-body forces and the comparison with experimental data. On the other hand, 
reactions induced by weak probes, like neutrino interactions with nucleonic matter, are relevant to astrophysics 
and can be calculated with few-body techniques. In this case, since often no experiment is possible, ab-initio 
predictions provide valuable input for astrophysical simulations. 



1 Introduction 

Excitations of nuclei with external electroweak probes are 
a fundamental tool to get information about nuclear dy- 
namics. If the probe is energetic enough, it can excite the 
system such that disintegration channels are open, and nu- 
cleons or compound fragments are emitted. In this case one 
would talk about an electroweak induced reaction. This is 
almost always the case if one considers light nuclei, as they 
possess no or very few bound excited states. In general, 
since a reaction implies states belonging to the continuum 
spectrum of the nucleus, exact direct calculations are dif- 
ficult and become prohibitive with increasing mass num- 
ber A. One is then forced to introduce some approxima- 
tions, that can be possibly validated against experiment. 
Such an approach is however not always possible. The first 
reason is that in nuclear physics the way the nucleons in- 
teract with each other, described by the Hamiltonian, is the 
object under investigation. Thus, comparing an approxi- 
mate calculation with experiments can pollute a clear in- 
terpretation on the reliability of the Hamiltonian used. A 
second reason is that often reactions that are relevant to 
astrophysics cannot be tested experimentally. That is why 
exact ab-initio calculations are needed. With this terminol- 
ogy we denote calculations where all nucleons are treated 
as degrees of freedom (ab-initio approach) and where the 
quantum many-body problem is solved numerically with- 
out introducing approximations (exact). Such calculations 
can be done for light nuclei. They can provide a valuable 
test on our present knowledge of nuclear forces via a com- 
parison with experimental data, and can also give useful 
predictions to be used e.g. in astrophysical simulations, 
when experiments are not possible. 

The interactions among nucleons are governed by quan- 
tum cromodynamics (QCD). In the low energy regime rel- 
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evant to nuclear structure and reactions, QCD is not per- 
turbative and thus difficult to solve. Nuclear theory has 
come up with a series of models to describe the interac- 
tion among nucleons in terms of effective degrees of free- 
dom, protons and neutrons. Several approaches have been 
investigated, starting from purely phenomenological ones, 
or using meson exchange theories and more recently via 
the systematic approach of effective field theory ||1|2| . In 
fact various realistic two-body potentials are available, that 
fit N-N scattering data with high accuracy. On the other 
hand, since the nuclear potential has clearly an effective na- 
ture, it is in principle a many-body operator. Furthermore, 
it is well known that most realistic potentials do not ex- 
plain the binding energy of three-body nuclei, unless three- 
body forces (3NF) are introduced. Today, aided by effec- 
tive field theory, where three-nucleon forces arise naturally 
and consistently with two-nucleon forces, a new debate is 
taking place regarding the role of 3NFs and how to con- 
strain them. For the determination of a realistic three-body 
potential or to discriminate among different models one 
needs to find A> 3 observables that are sensitive to 3NFs. 
An important activity in this direction has taken place in re- 
cent years, with accurate calculations of bound-state prop- 
erties of nuclei of increasing mass number A II3I4L A com- 
plementary approach should be direct to the study of light 
nuclei reactions. We have decided to draw our attention 
towards this second research line, concentrating on elec- 
troweak induced reactions. 



The paper is structured in the following way. We will 
first present the main calculational techniques used to tackle 
the problem of electroweak reactions with light nuclei in 
Sec. 121 Then, we will separate the discussion into electro- 
magnetic observables in Sec. [3] and weak observables in 
Sec.|4] Finally, we will conclude in Sec. |5] and present an 
outlook for the future. 
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2 Calculational Techniques 

In the investigation of electroweak reactions one funda- 
mental ingredient is the nuclear response function to the 
excitation of an external probe. In an inclusive process it is 
defined as 




where w and q are the energy and momentum transferred, 
I !?o/ /) and £0// denote initial and final state wave functions 
and energies, and finally O represents the excitation opera- 
tor, which is different according to the external probe. The 
5-function ensures energy conservation. From Eq. ([1]) it is 
evident that in principle one needs the knowledge of all 
possible final states excited by the electromagnetic probe, 
including of course states in the continuum. Thus, in a 
straightforward evaluation one would have to calculate both 
bound and continuum states. The latter constitute the ma- 
jor obstacle for a many-body system, since the full many- 
body scattering wave functions are not yet accessible for 
A > 3. A way to circumvent the problem is to use the 
Lorentz Integral Transform (LIT) method [5 '6] . Within the 
LIT one considers instead of Rgito, q) an integral transform 
JLoicr, q) with a Lorentzian kernel defined for a complex 
parameter cr - ctk + i ctj by 

£oicr, q) = r da. "'^""f ^ = {W^J^^J . (2) 

The parameter cr/ represents the resolution of the trans- 
form and is chosen to have finite value (ctj + 0). The basic 
idea of considering £,0 lies in the fact that it can be evalu- 
ated from the norm of a function ^¥^^^, which is the unique 
solution of the inhomogeneous equation 

(i/-£o-cr)|^^,) = 0(^)1 ^'o). (3) 

Here H is the nuclear Hamiltqnian. The existence of the in- 
tegral in Eq. (|2| implies that has asymptotic boundary 
conditions similar to a bound state. Thus, one can apply 
bound-state techniques for its solution. 

In all the results presented here, if not indicated differ- 
ently, in order to solve Eq. (O we expand ^'o and ^'^^ in 
terms of antisymmetrized hyperspherical harmonics (HH) 
121. The expansion is performed up to a maximal value of 
the grand-angular momentum K^ax- We improve the con- 
vergence using an effective interaction defined within the 
hyperspherical harmonics (EIHH) |8|, where the starting 
potential is replaced by an effective one constructed via 
the Lee-Suzuki method. 

The response function Rq is then obtained by numer- 
ically inverting the integral transform in Eq. For the in- 
version of the LIT various methods have been devised |9 10]. 
In particular the concept and methods of the inversion have 
been recently clarified and discussed extensively in Ref. IfTTll . 
in response to some criticisms. Here, we will take care in 
pointing out the error bands for the calculated theoretical 
curves due also to the numerical inversion procedure. 

In the next sections we will present recent applications 
of these calculational techniques to electromagnetic and 
weak reactions with light nuclei. 



3 Electromagnetic Observables 

3.1 Photo-absorption Reactions 

In recent years there has been an increased interest in the 
photo-disintegration of light nuclei, both in theory and ex- 
periment. From the theoretical side, a lot of effort has been 
addressed to the investigation of such process with the LIT. 
Several light nuclei have been studied with mass number 
up to A — 1 1121131 . This has stimulated new experimen- 
tal activity in Lund ifTSll . for example. The most system- 
atic calculations with the LIT have been performed on the 
'^He target, where a variety of realistic potentials have been 
tested (see e.g. [14], and reference therein). 

One important reason for concentrating on "^He is that 
it is the lightest nucleus that already possesses features of 
heavier systems (e.g. large binding energy per nucleon). 
Therefore "^He can be considered as an important link be- 
tween the classical few-body systems and more complex 
nuclei. Furthermore, continuum properties of "*He are an 
ideal testing ground for microscopic three-nucleon forces, 
which are typically fitted to bound state properties of few- 
body systems. 




Fig. 1. Feynman diagram of the photo absorption process. 

The total photo-absorption cross section, represented 
by the Feynman diagram in Figure [1] is given by 

(TyioS) = \Ti^aujRoX(jS) , (4) 

where a is the fine structure constant and Ri,_ is the re- 
sponse function in the unretarded dipole approximation 

RuXoi)^ ^ dWi\(^'j\D,\¥j\5(Ef-EQ-oS), (5) 

with 

being the unretarded dipole operator. Here, and zy^ repre- 
sent the third component of the isospin operator and of the 
spatial coordinate of the ^-th particle in the center of mass 
frame. This is a typical example of an electromagnetic in- 
duced reaction that can be studied with the LIT. 

Photo-nuclear processes can be of great importance, 
since exchange currents, which are connected to the nu- 
clear interaction by current conservation, play a decisive 
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role in the cross section. Within the unretarded dipole ap- 
proximation, the dominant part of the exchange current is 
taken into account implicitly via the Siegert theorem, par- 
ticularly for photon energies oj below 50 MeV. This has 
been shown for the deuteron |16| and is expected to be 
valid for '^He as well. 
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Fig. 2. (a) Total "^He photo-absorption cross section CTy for AV18 
(dark blue) and the^EFT N^LO (light blue) potentials; (b) cTj in- 
cluding the 3NFs UIX (red) and the ^EFT N^LO (pink) in com- 
parison with with some experimental data: box from P211, squares 
from L22J . and circles from L23J . 



The first realistic calculation of the He photo-absorption 
was performed with the LIT and EIHH method 1 17] using 
the Argonne Vm (AY 18) 1.18 1 potential augmented by an 
Urbana IX (UIX) lfT9l three-nucleon force. Later on, a sec- 
ond realistic calculation was performed with the LIT and 
the no-core shell model (NCSM) |20| using chiral two- 
and three-body forces IT]. It is very interesting to com- 
pare these two calculations in the following way. First one 
can look at the predictions from two-body forces only. Fig- 
ure |2a) shows in fact the result for o-y(aj) with the AVI 8 



and the ;i^EFT force at N^LO |[T|. One observes that, even 
though both forces reproduce very well the nucleon-nucleon 
(NN) experimental phase shifts (with a pei" datum of 
about 1), the prediction of this many -body observable is 
different, due to the different off-shell properties of the 
two-potentials, which are not fixed by the NN phase-shifts. 
Figure |2b) shows, instead, cry{oj) calculated from Hamil- 
tonians that include a three nucleon force: in one case the 
UIX and in the other the^EFT force at N^LO. We observe 
that the two predictions are very similar, even though the 
3NF have been constrained differently (the UIX on the tri- 
ton binding energy and on saturation properties of nuclear 
matter and the chiral N^LO on s- and p-shell nuclei [4]). 

One should note that the theoretical error bands shown 
in Fig.|2]take into account the numerical error of the LIT 
inversion. Furthermore, the AVI 8-1- UIX curve has a larger 
width due to the fact that an extrapolation was done to an 
infinite model space for the LIT itself ilTl . 

It is particularly interesting to observe the relevance of 
3NFs. Due to 3NFs one observes a reduction of the peak 
height by less than 10% both for the more phenomenologi- 
cal potentials as well as in the chiral approach. One should 
note, though, that larger 3NF effects were found with the 
UIX in [17] at higher energies than those shown in Fig.|2] 

In Fig. |2b) we also compare to some selected exper- 
imental data (for more data see [17] and [20]): the box is 
from |21 1, where the peak cross section is determined from 
Compton scattering via dispersion relations; squares from 
[i22| and dots from [23] show measurements of the cross 
section in a broader energy range. It is evident that the var- 
ious experimental cross sections are quite different exhibit- 
ing maximal deviations of about a factor of two. While the- 
ory is relatively accurate, considering the sum of the two 
theoretical bands in Fig.|2b) as an estimate of the total the- 
oretical error, the experimental situation is unsatisfactory 
and does not allow one to draw any definitive conclusions 
on the importance of 3NFs to describe this observable. 



3.2 Electron Scattering Reactions 

The information about the nuclear dynamics that one can 
get from electron scattering processes is complementary to 
the one obtained with real photons and much richer. While 
the energy of the photon fixes the momentum transferred 
to the target, an electron can transfer energy o) and mo- 
mentum q independently. As q - |q| increases one can 
probe properties of the correlations among nucleons from 
the long- to short-range. 

Various inclusive experiments have already been per- 
formed in the past on light nuclei. Due to the low atomic 
number it is possible to study the longitudinal and the trans- 
verse responses without the ambiguities created by Coulomb 
distortions affecting heavier systems. 

We report on recent studies of the the inclusive inelastic 
electron scattering reactions on light nuclei with A -A and 
3 on a broad energy range obtained using the LIT method. 
This is, in fact, the only way to allow investigations beyond 
the three- and four-body (for A - A) break up thresholds. 
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Fig. 3. Feynman diagram of the electron scattering process at the 
one-photon exchange level. 



In the one-photon-exchange approximation (see FiglSj, 
the inclusive cross section for electron scattering off a nu- 
cleus is given in terms of two response functions. 



dQda> 



- CTm 



91 



—Rdo), ^) + 1 ^ + tan - 1 RAoi, q) 



(7) 

where ctm denotes the Mott cross section, Q- - -q^ - q^- 
Ci?, and the electron scattering angle. The longitudinal 
and transverse response functions, Rl{cj, q) and Rt{cl>, q), 
are determined by the transition matrix elements of the 
Fourier transforms of the charge and the transverse cur- 
rent density operators, respectively. Experimentally R^ and 
Rt can be disentangled by using a Rosenbluth separation 
method. 

We will start presenting a study of the longitudinal re- 
sponse function for "'He and and then continue with the 
transverse response function for ''He and the A - 3 nuclei. 
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3.2.1 The Longitudinal Response Function 
The longitudinal response function is given by 

R,(aj,q) = ^ \{'Pf\p(qWo)f disf + ^ - Eo - 

^ (8) 
where M is the target mass. The charge density operator p 
is defined as 



P(q) ^^^k exp [iq ■ rt] , 



(9) 



where is the nucleon charge in unit of 43 and r^; is the 
coordinate of the k-th particle in the center of mass frame. 
In order to take the finite size of the nucleon into account, 
form factors are introduced, such that the nucleon charge 
takes the following form 



ek 



— 7^G^(e') + — r^G^e^) , (10) 
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Fig. 4. Longitudinal response function for ''He at various q: 
PWIA (dotted) versus the calculation including the full FSI 
(solid) from the AV18+UIX potential. Data shown are from 
Bates L30J (squares), Saclay L31J (circles). 



where G^^" are the proton and neutron electric form fac- 
tors. We use the usual dipole fit for them ll24l . 



e is here factorized in ctm- 



Since the longitudinal response Rl is much less sensi- 
tive to meson exchange effects than the transverse response 
Rt, the use of a simple one-body density operator allows 
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one to concentrate on the nuclear dynamics generated by 
the potential. In fact, for low q two-body operators in 
are only of fourth order in effective field theory counting 
(N^LO) II25I . and their contribution is negligible up to ^ 
300 MeV/c. In recent works |26 27 1 we have studied the 
longitudinal response function Ri((L),q) for ^He at con- 
stant momentum transfers q < 500 MeV/c. We consider 
the non-relativistic charge density operator of Eq. (|9]l, and 
decomposed it into isoscalar and isovector Coulomb mul- 
tipoles fTT\. Each multipole is then used as source in the 
right-hand-side of Eq. (O which is solved with the EIHH 
method. After inversion of each multipoles and a subse- 
quent sum of all of them, one obtains results that can be 
directly compared with experiment. The accuracy of the 
convergence for the LIT is found to be at a percentage level 
for all the results shown in the next paragraphs, including 
the numerical error of the inversion. 

It is first interesting to investigate the role of the final 
state interaction (FSI) in the electron scattering reaction. In 
ll28l it was found that the FSI is very important for Rl in 
the quasi-elastic peak of "^He. To investigate that we com- 
pare our LIT results using the AV 1 8+UIX potential both in 
the initial and final state with a plane wave impulse approx- 
imation (PWIA) calculation. The PWIA result is obtained 
under the hypothesis of one outgoing free proton with mass 
m and a spectator (A-l)-system with mass M^: 

«rv.,)=/.p,.(p).(.-*i2^-^-.). 

(11) 

Here n(p) represents the proton momentum distribution 
and e the proton separation energy. To obtain a PWIA curve 
we have used the n(p) of 1291 . In Fig. |4] the results of 
RL{(L),q) with the AVISh-UIX potential at various q are 
shown and compared to data. In all cases one finds that 
the FSI effects are very large and essential for reaching 
agreement with experiment, confirming what was previ- 
ously found. One has to note, that the PWIA fails particu- 
larly in the q.e. peak and at low oj, while it is better on the 
tail. 

The second fact one can investigate is the role of three- 
nucleon forces on this observable. To this aim, we can for 
example compare calculations with a two-body force only, 
like the AV18 potential, with the one with AVlS-hUIX. In 
Fig. |5j we show this comparison for some values of q to- 
gether with existing experimental data. One sees that the 
3NF results are closer to experiment, this is particularly 
evident at ^ = 300 and 350 MeV/c. The 3NF effect is gen- 
erally not very large, but always leads to a quenching of the 
strength. For some q values the 3NF quenching is compa- 
rable to the size of the experimental error bars, particularly 
for the data from Ref. |30|. The largest discrepancies with 
data are found at ^ = 500 MeV/c. While the height of the 
peak is well reproduced by the result with 3NF, the width 
of the experimental peak seems to be somewhat narrower 
than the theoretical one. On the other hand one has to be 
aware that relativistic effects, which we are not including, 
are not completely negligible at such momentum transfers 
(we will comment on that later). 



It is worth looking at the 3NF effect when even lower 
momentum transfers are considered, even though for some 
kinematics there are no experimental data to compare to. 
In Fig. |6] we show Ri at three momentum transfer 50 < 
q < 200 MeV/c. We find a large quenching effect due 
to 3NFs, which is strongest at lower q. We investigate it 
with two different 3NF models: the UIX and the Tucson 
Melbourne (TM') [32J, when starting from the same two- 
nucleon AVI 8 potential. While the UIX force contains a 
two-pion exchange and a short range phenomenological 
term, the TM' force includes two-pion exchange terms where 
the coupling constants are taken from pion-nucleon scat- 
tering data consistent with chiral symmetry. The cutoff of 
the TM' force has been calibrated on triton binding energy, 
when used in conjunction with the AVI 8 NN force. Thus, 
the TM' is not adjusted also to nuclear matter properties as 
the UIX. With a cutoff A - Amm„, where m„ is the pion 
mass, we obtain the following binding energies 8.47 MeV 
(^H) and 28.46 MeV (^He) for the AV18-HTM' potential. 
We would like to emphasize that the '*He binding energy is 
practically the same as for the AVI8-1-UIX case (as already 
found in |33 1). 

Figure|6]shows that the increase of 3NF effects with de- 
creasing q is found for both the UIX and TM' force model. 
The quenching of the strength is rather strong, reaching 
40% for ^ = 50 MeV/c. Moreover it becomes evident 
that the difference between the results obtained with the 
two 3NF models increases with decreasing q. One actu- 
ally finds that the shift of the peak to higher energies in 
the case of UIX generates for Ri a difference up to about 
10% on the low energy sides of the peaks. This is very in- 
teresting, since it represents the case of an electromagnetic 
observable whose prediction depends on the choice of the 
3NF. In light of these results it would be very interesting 
to repeat the calculation with EFT two-and three-body po- 
tentials 1 1 2 1 . While to our knowledge there are not avail- 
able measurements at the low momentum kinematics like 
^ = 50 and 100 MeV/c, recently some data have been taken 
at ^ ^ 200 MeV/c |34|. They are shown in Fig|6] While 
one finds a satisfactory agreement between the AV18-i-UIX 
result and data beyond the peak, one observes a non negli- 
gible discrepancy at the peak itself. Moreover the error bars 
are quite large. It would be highly desirable to have precise 
measurements of Ri at low q. This could serve either to fix 
the low-energy constants (LEC) of the effective field the- 
ory three-body force or to possibly discriminate between 
different nuclear force models. Here, it is worth mention- 
ing that our calculations have stimulated new experimental 
activity in Mainz ll35l . More experiments on the same ob- 
servable also planned at Jefferson Laboratory (E05.110 at 
Hall A) at large momentum transfer 

It is interesting to compare the '*He longitudinal re- 
sponse at low momentum transfer with the corresponding 
calculations done for A - 3 nuclei published in lf36l . Also 
for A = 3 nuclei a reduction of the peak due to the 3NF 
was observed and it was found to be more pronounced at 
low momentum transfer. Though the percentage effect is 
different for A = 4 and A = 3 nuclei. In case of the a- 
particle the 3NF force is quenching the quasi-elastic peak 
of about 20% for q = 200 MeV/c. For ^H, instead, the 
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quenching effect of 3NF was found to be of about 10% 
only at ^ = 174 MeV/c at energies between 10 and 20 
MeV. This can be seen in Fig|7] where theoretical results 
with the same potentials used for "^He are shown. In addi- 
tion, a calculation with the configuration space BonnA 1371 
potential augmented by the TM' force was done. Compar- 
ing the three cases in FigQwhere a 3NF was included, one 
finds that the result is very stable. 

Summarizing, in case of "^He two features are evident 
with respect to ^H: (i) there is a larger eff'ect of 3NFs and 
(ii) the two diff'erent three-body force models lead to two 
diff'erent curves. This shows that '^He is a more interesting 
target to study 3NFs, since their eff'ect is amplified. This 
is expected since it is a denser systems and the ratio be- 
tween the number of triplets and of pairs goes like (A-2)/3, 
therefore doubles from ^'H to '*He. 

Having the full response function at our disposal, we 
can also calculate some integral properties to compare with 
other calculations found in literature. The classical sum 
rule of Rl{(^, q), that has been much discussed in the lit- 
erature, is the Coulomb sum rule (CSR) t38 |- It consists in 



the integral of the inelastic longitudinal response function 
CSR(,)= [L-^p^ 

= Z + Z(Z - l)f„„(q) - Z^\F(q)f, (12) 

that can be connected to the number of protons Z and to 
fppiq)^ which is the Fourier transform of Ppp{s), i.e the 
probability to find two protons at a distance s. Here F{q) is 
the nuclear elastic form factor. 

We can first look at the 3NF eff'ects on the CSR. In 
Table [1] we show the CSR for different momentum trans- 
fers obtained with the AV18 and the AVlS-hUlX potential. 
Clearly, the UIX force leads to a quenching of the CSR, 
especially at low momentum transfer, consistent with what 
was observed in the response function. At ^ = 50 MeV/c 
the quenching on the overall strength is of about 20%. 

We can now concentrate on fppiq), since it contains the 
interesting physical information about the proton-proton 
correlation function and since we can compare it with pre- 
vious calculations. In Fig. [8] fppiq) is shown in compar- 
ison with the results of Ref. |39| obtained with the same 
potential. In Fig. [8] it is interesting to observe the perfect 
agreement of the two results at low q, where the one-body 
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Fig. 6. Longitudinal response function for low-momentum trans- 
fer with the AV18 two-nucleon force only (dashed), and with the 
addition of the UIX (solid) or the TM' (dashed-dotted) three- 
nucleon force. 



charge is used. We can also see that the contributions of 
two-body terms, called meson exchange currents (MEC), 
to the charge density operator become negligible below 
q = 300 MeV. 

Besides the non relativistic expression p(q) = pm((l) 
of Eq. (|9]i, the first contribution due to relativistic effects 
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Fig. 7. Longitudinal response function of for q = 174 MeV/c 
with the AVI 8 two-nucleon force only (dashed), and with the ad- 
dition of the UIX (solid) or the TM' (dashed-dotted) and with 
the Bonn two-body potential with the TM' 3NF (dotted). Calcu- 
lations are taken from [36 1. 

Table 1. CSR of Eq. l ll2t for different momentum transfers ob- 
tained with the AVI 8 and the AV18+UIX potential. 
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body density operator (dots), and with the inclusion of MECs 
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lead to a one-body operator given by 

^2 



P«c(q) = > , ek exp [/q ■ r^] + 



(13) 



Am? 



cT/i ■ (q X Pi) exp [/q ■ Ti] . 



where 



1 +T 



(14) 



with G^" being the proton and neutron magnetic form fac- 
tors. Here crk and p^^ are the spin and the momentum of 
the A:-th particle, respectively, and m is the mass of the 
nucleon. The first term in Eq.(fT3ll is denominated by the 
Darwin-Foldy term, while the second is a spin-orbit one- 
body term. In all previous Figures for '^He we have ne- 
glected relativistic contributions (RC). 
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Fig. 9. Longitudinal response of '^He with the non relativistic 
(NR) one-body current of Eq. (|9} and with the Darwin-Foldy rel- 
ativistic correction (RC) (see text), for q = 450 and 500 MeV/c. 



the Darwin-Foldy one body relativistic correction in Fig|9] 
where its effect at ^ = 450 and 500 is shown. It leads to 
a decrease of the cross section in the quasi-elastic peak by 
5 or 6%. The spin-orbit correction was found to be small 
in the case of A = 3 |36|, and is probably small in "^He as 
well. Other relativistic effects, like wave function boosts 
have been found to be important in Ref. |36|. We plan to 
investigate their effect on "^He in the near future. 



3.2.2 The Transverse Response Function 
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Fig. 10. Transverse response function Rj of ''He slI q = 100 and 
q = 300 MeV/c calculated with the MT potential: one-body cur- 
rent (dashed curve) response and one- plus two-body current re- 
sponse (solid curve) in comparison with experimental data from 
Bates 1^ and Saclay (311. 



The investigation of the transverse response function 
in the electron scattering reaction is more involved, as it 
requires the knowledge of the transverse current operator. 
The transverse response function is defined as 



While the inclusion of two-body currents in a realistic 
calculation of Rl requires some labor, it is easier to in- 
clude first order relativistic corrections. We have included 



RT(A>,q) 



^\{^f\jT(q)mf^(Ef + ^-Eo 



(15) 
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The transverse electromagnetic current density of a nu- 
cleus can be decomposed in a superposition of one- and 
many-body operators 

jr(q) =j(i)(q)+j(2)(q) + --- ■ (16) 

The non-relativistic one-body current density consists of 
convection and spin currents 

j(i)(q)=j(i)(q)+j(i)(q), (17) 

with 

j(i)(q) = ^^ei{Pi,exp[;q-ri]}, (18) 

k 

j(i)(q) = ^ ((Tk X q) exp [iq ■ r^] , 

k 

where and fik are the same as in Eq. ( fTOl ) and (fT4l i. re- 
spectively. Even in a non-relativistic approach, the two- 
body term j(2) does not vanish. In fact, a momentum and/or 
isospin dependent two-body interaction V requires the ex- 
istence of a two-body current operator, in order to satisfy 
the continuity equation, which in coordinate space, reads 

V-j(2)(x) = -/[y,p(x)], (19) 

where p is the Fourier transform of the one-body charge 
in Eq. (|9]l- As is well known, this relation is not sufficient 
to determine uniquely the two-body current. Therefore one 
needs in principle a dynamic model for the nuclear poten- 
tial which reveals the underlying interaction mechanism. 
Such a model is supplied, for example, by the meson ex- 
change picture of the NN force. But even for a phenomeno- 
logical potential a consistent MFC can be constructed by 
expanding the potential in meson exchange-like terms. 

Investigations of Rt with the LIT have been done for 
A = 4 and A - 3 nuclear targets. For the a-particle only 
calculations with semi-realistic forces have been done BOl 
so far, using central potentials with spin-isospin depen- 
dence, like the Malfliet-Tjon (MT) |4U. In Figure [10] we 
show our main results for q = 100 and q - 300 MeV/c. 
We present a calculation where only the one-body cur- 
rent operator is used, called the impulse approximation 
(lA), in comparison to one where consistent MECs are 
included. Quite strong effects of the MFC are found at 
low-momentum transfer, where at energies of about 100 
MeV, the two-body currents lead to an enhancement of the 
strength by a factor of 4 with respect to the lA. We also ob- 
served that at this low momentum transfer the spin current 
in Eq. (fTST l does not dominate strongly: it only leads to half 
of the total lA strength |40|. That is why Rt at low q shows 
sensitivity to the convection current in Eq.dTSll as well as 
to the MFC. Unfortunately, to the best of our knowledge, 
no experimental data are available for these kinematics. At 
larger momentum transfer, where data are found in liter- 
ature, only a small effect due to the MFC is found. Here 
we observe that the spin current strongly dominates the 
strength at these kinematics. Moreover, the weak MFC ef- 
fect is possibly due to the semi-realistic nature of both the 



potential and the consistent MEC, which included only fic- 
titious scalar meson exchanges and no n- and p-exchange 
currents. 

Realistic calculations with the LIT have been done, in- 
stead, on the A - 3 body systems (■'H and ^He) [421 and 
P3l . The interest of this kind of study is not only to inves- 
tigate the role of 3NFs. The electron scattering reaction on 
the ^He nucleus, in fact, is of particular importance, since 
it can also serve as effective neutron target. 

In (47] the BonnA [37] two-body potential was chosen, 
due to its meson-exchange character. In this way, meson 
exchange currents could be determined relatively uniquely. 
The Hamiltonian was also augmented by an additional three- 
nucleon force, the TM'. MEC of n- and p-exchange na- 
ture have been used. In a subsequent calculation |43|, the 
AV18H-UIX force model was used, instead. In this case 
consistent MEC for the AVI 8 potential were introduced 
using the Arenhovel and Schwamb model [44], where an 
expansion of the potential in terms of fictitious meson ex- 
change is done. Theoretical results were obtained for dif- 
ferent current operators. A one-body current as in Eq. ( fTsT i 
was used with the addition of one-body relativistic terms 
[43 [ and MECs. In Fig. [TT] we compare theoretical results 
from the AVISh-UIX potential with the experimental data 
from Retzlaff et al. Il45l . Both the experimental data and 
the theoretical curves are plotted as a function of the in- 
trinsic energy u;,-,,,, i.e. the recoil energy ^ written in the 
formula (fTSl l has not been added. In the case of ''He one 
has a rather good agreement of theoretical and experimen- 
tal results at the higher q values. The MEC contribution 
is essential to reach agreement with the experiment and is 
pretty large at the threshold. At q = 174 MeV/c the theo- 
retical prediction for Rj underestimates the data below 1 1 
MeV. In the case of ^H the comparison of theory with re- 
spect to experiment looks worse, with a systematic under- 
estimation of data, which gets larger at the lowest ^-value. 
Relativistic contributions are small, but visible already at 
q = 323 MeV/c. In [43 [ it was also observed that different 
force models, like the BonnA+TM' and the AV18-I-UIX 
lead to very similar result. 



4 Weak Observables 

In this section we deal with reactions induced by neutrino 
interactions on nucleonic matter Such processes are rele- 
vant to astrophysics, for example in the context of core- 
collapse supernovae. The present theory of supernovae ex- 
plosion still presents open questions about the explosion 
mechanism itself and also regarding late nucleosynthesis 
processes. To better understand these aspects, more micro- 
scopic calculations and improved simulations are needed. 
Since, often no experiments are possible to simulate the 
conditions of a supernovae, few-body calculations can pro- 
vide valuable insight on many relevant reactions, that can 
then be used as input for hydrodynamic simulations of the 
explosion. We will first discuss about inelastic neutrino re- 
actions of "^He and secondly about neutrino-antineutrino 
bremsstrahlung in neutron matter. 
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Fig. 11. Comparison of experimental 1451 and theoretical results for the threshold ^^He (upper panels) and (lower panels) transverse 
response function Rj with the AV18+UIX. Theoretical curves show different current operators: relativistic one-body current (dashed), 
relativistic one-body current + MEC (solid), non relativistic one-body current + MEC (dotted). 



4.1 Neutrino Scattering 

Core-collapse supernovae are believed to be neutrino driven 
explosions of massive stars. At the end of its life a massive 
star (M > 8 Mq, where Mq is the solar mass) presents 
a typical onion-like structure with layers containing ashes 
of different burning stages, ending up with an iron core 
at the center. As the iron core grows and becomes grav- 
itationally unstable, it collapses until nuclear forces halt 
the collapse and an outgoing shock wave is created. The 
most likely scenario is that the shock wave that propagates 
in the outer layers is stopped at some point due to energy 
loss via the dissociation of nuclei into mainly free nucleons 
and o'-particles. One mechanism that can trigger a revival 
of the shock is given by interactions of neutrinos coming 
out of the proto-neutron star (the remnant at the center of 
the iron core) with matter in the outer layers. In such a de- 
layed explosion scenario, the neutrino energy transfer oc- 
curs mainly via elastic scattering of vs with matter (free nu- 
cleons and nuclei). Nevertheless, other reactions can take 
place and contribute to the energy transfer. One of them is 
the inelastic excitation of light nuclei, like ^H, ^He, ^H, and 
"^He f49l, fSO]. The hot dilute gas above the proto-neutron 
star and the region below the accretion shock contains up 
to 70% of "*He. That is why the inelastic excitation of "^He 



via neutrino scattering is an important process to study in 
such a context. 
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Fig. 12. Feynman diagram of the inelastic neutrino scattering off 
a nucleus, where all neutrino flavors can be involved. 



Neutrinos are produced in flavor equilibrium inside the 
proto-neutron star The characteristic temperatures of the 
emitted neutrinos are though different, due to the differ- 
ent reactions they undergo. They are about 6-10 MeV for 
Vfi^riVf^ j), 5-8 MeV for Vg and 3-5 MeV for v^,. Thus, mainly 
Vf^^T and v^ j can carry more than 20 MeV of energy, which 
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is the amount needed to dissociate "^He via an inelastic neu- 
tral current interaction. The Vg and Vg can also contribute to 
a dissociation of the a-particle via a charged current inter- 
action, but in a less important way due to their lower en- 
ergy. Even though a number of studies have been done in 
the past to estimate these cross sections, the first ab-initio 
few-body study was done by Gazit and Barnea |46| using 
the LIT and EIHH methods. A calculation was done for 
the energy dependent inclusive inelastic cross section for 
4He(v^, v'Jpi, ^He(v^, v'Jpi, ^ne(v„ e+)% and ^He(ve, 
e")^X, where x = e,fi,T and is the final state of A- 
nucleons with Z-protons. The corresponding Feynman di- 
agram is represented in Fig. [12] 
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Fig. 13. Axial quadrupole transition in the neutral neutrino scat- 
tering off ■*He as a function of the excitation energy a>. Calcula- 
tions are performed with the AVI 8 and the AV18+UIX potential 
and are taken from Ref. [47]. 



Table 2. Temperature averaged neutral current inclusive inelas- 
tic cross section per nucleon in case of ^He as a function of the 
neutrino temperature. Data in [10"*^ cm-] taken from |46|. 
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The calculation was done using the AVI 8 two-body 
force augmented by the UIX three nucleon forces. The 
effect of the 3NF on astrophy sic ally relevant observables 
is an interesting one to investigate. In Figure [13] the re- 
sponse function of the dominant multipole (the axial E2) 
in Eq.(l20b for the neutral interaction cross section is shown 
as a function of the transferred energy co for the AVI 8 and 
AV18+U1X potential. The numerical accuracy of these cal- 
culations is of the order of 1%. It is interesting to note that 
also in this case one observes a quenching effect due to the 
UIX force. This is analogous to what was found for the 
electron scattering observable, even though we are dealing 
with two different cuiTent operators. 

The effect of MEC was also investigated using an hy- 
brid approach: the weak current was constructed within 
EFT and then used with wave functions obtained from the 
traditional potentials AVI 8 and UIX. The axial current in- 
troduces a dependence of the MEC on a new low-energy 
constant, denoted with d,-, which was calibrated on the tri- 
ton half life. The effect of the 3NF and the MEC on the 
temperature averaged neutral cross section 



After a multipole decomposition of the weak current, 
the cross section for the neutrino scattering process can be 
written as [46 1 

-klF"(Zf,kf) (20) 



Ida" 



,«'^//v(v) 27, + 

f CO 

X i 2 [XcRcj + XlRlj - XcL^^Rc'^L'] 



q 

sin 6 (5 e - w H 

;o I 2M 
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where G is the Fermi weak coupling constant, /, is the ini- 
tial angular momentum of the nucleus, kj is the momen- 
tum of the outgoing lepton, Zj is the charge of the residual 
system and F is the Coulomb function. The superscript a 
denotes the isospin component, with a - Q for the neutral 
current and a - +\ for the charged current. Xo^o2 is the 
lepton kinematical function BSll and RoiOi is the response 
function to the operators 0\ and 0\, where the notation 
Cj,Lj, Mj and Ej stands for Coulomb, longitudinal, mag- 
netic and electric multipoles of order J. 



(21) 



is shown in Table [2] Here the (. . . )r symbol indicates the 
temperature average 



r°° r dcr° 

doj dkif(T,,kd—, (22) 
Jm J dkf 



which is obtained by integrating the cross section in energy 
and assuming that neutrino spectra are described by Fermi- 
Dirac distributions / with temperature Ty. Here fc,- is the 
initial momentum of the neutrino. One observes that the 
low-energy cross section is very sensitive to the nuclear 
force model, while MFCs are contributing in a negligible 
way. This is explained by the fact that the multipole where 
MFCs are expected to contribute the most is the Gamow- 
Teller, which is suppressed in ^He. 

One could ask the question whether the inelastic exci- 
tation of nuclei other than the ff-particle can contribute to 
the neutrino energy transfer to the shock wave in super- 
novae. In order to tackle that question, one not only needs 
to know the inelastic cross section, but also the composi- 
tion of the nucleonic matter. The equations of state (EOS) 
typically used in supernovae simulations include neutrons. 
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Fig. 14. Neutrino energy loss dEy/dx for inelastic excitations of A = 3, 4 nuclei as a function of the matter temperature T and at density 
of 10'^ g/cm^. Neutrino energy spectra are assumed to be Fermi-Dirac distributions with temperature Ty = 6 MeV. Different proton 
fractions are shown. Figure taken from Ref. L49J . 



protons, a-particles and some representative heavy nuclei. 
In a recent study by O'Connor et al. 149 i the virial EOS 
was extended to include also and ^He nuclei and pre- 
dict A = 3 mass fractions near the neutrino sphere in su- 
pernovae. Neutrino cross sections of A = 3 inelastic exci- 
tation were calculated with the LIT and EIHH, using the 
same potentials and MEC as done for the a-particle. An 
important fact was found: even though the a-particles are 
more abundant, the loosely bound A = 3 nuclei dominate 
the energy transfer at low-density via inelastic reactions. 
MECs were also found to be important: they contribute up 
to 16% in the thermal averaged cross section at tempera- 
ture of 1 MeV. This is explained by the fact that the thresh- 
old cross section (which is dominant at low Ty) is largely 
influenced by the the axial dipole transition, where MECs 
are important. 

In Figure [T4l the neutrino energy losses. 



dEy / 0\ 

i=3H ^He.-'He 



(23) 



are shown. Here Ey is the neutrino energy, nh is the baryon 
density, x; - ^ is the mass fraction and (wcr),- j-^ is the en- 
ergy weighted temperature averaged cross section (same 
as in Eq. (l22l l with an energy weight). The neutrino en- 
ergy loss for inelastic excitations of A = 3,4 nuclei was 
calculated as a function of the matter temperature T ~ A 
MeV and at density of about 10'^ g/cm-', which are typ- 
ical conditions near the neutrino-sphere. Neutrino energy 
spectra are assumed to be Fermi-Dirac distributions with 
temperature Ty - 6 MeV. Different proton fractions Yp - 
0.1, 0.3, 0.5 were investigated. One observes that A - 3 nu- 
clei contribute significantly to the neutrino energy loss due 



to inelastic excitations at T > 4 MeV. As such, they should 
be included in supernovae simulations. 



4.2 Neutrino-Antineutrino Bremsstrahlung 

Neutrino processes involving two nucleons play a key role 
in the physics of core-collapse supernovae and neutron stars. 
In particular, neutrino-pair bremsstrahlung and absorption, 
A^A^ <-> NNvv, are important processes for the production 
of muon and tau neutrinos lISDl. Due to the fact that neu- 
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Fig. 15. Feynman diagram of the neutrino-antineutrino 
bremsstrahlung process on two neutrons. The interaction 
between the two nucleons is here represented by a one-pion 
exchange term (pion denoted with the dashed line). 



trinos interact weakly, the rates for neutrino emission and 
absorption are basically determined by the dynamic re- 
sponse functions of strongly-interacting matter. Supernova 
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explosions are most sensitive to neutrino processes near 
the proto-neutron star, at subnuclear densities p ~ po/10, 
where po - 2.8 x 10'"* g/cnr' is the saturation density, and 
where matter is neutron rich. 

The bremsstrahlung rate was calculated from NN phase 
shifts in Ref. ||52|| . Here, we report on our recent calcula- 
tion of the neutrino bremsstrahlung in pure neutron mat- 
ter, where we used chiral effective field theory (EFT) [5 3 J . 
In Figure [15] a Feynman diagram for neutrino- antineu- 
trino bremsstrahlung is shown. Note that this is not the 
only Feymnam diagram involved, as the v v pair can be 
produced from the ingoing left or right neutron, or from 
the outgoing left neutron, as well. If one considers all the 
four possibility, as done explicitly in Il54ll . one can observe 
that only non-central contributions to strong interactions 
are relevant in such process. Thus, in a chiral EFT lan- 
guage, tensor forces from pion exchanges and spin-orbit 
forces, are essential for the two-nucleon response. 

In supernova and neutron star simulations, the stan- 
dard rates for bremsstrahlung and absorption are based on 
the one-pion exchange (OPE) approximation M54I55I for 
the interaction between the neutrons, which is represented 
by the dashed line in Fig. [15] This is a reasonable start- 
ing point, since it represents the long-range part of nu- 
clear forces, and it is the leading-order contribution in chi- 
ral EFT 1 561. However, sub-leading non-central contribu- 
tions are crucial for reproducing NN scattering data |56| 
and might be relevant in neutrino bremsstrahlung as well. 
We systematically go beyond the OPE approximation and 
include contributions up to next-to-next-to-next-to-leading 
order (N^^LO) in chiral EFT. We find that shorter-range 
non-central forces significantly reduce the neutrino rates 
for all relevant densities. 

We follow the approach to neutrino processes in nu- 
cleon matter developed in Ref. |[57l, which is based on 
Landau's theory of Fermi liquids. This is a good approx- 
imation since the energies a> and momenta q transferred 
to the neutrinos are small in comparison to the momenta 
of the neutrons. In addition, we consider degenerate con- 
ditions, where the temperature is small with respect to the 
Fermi energy. The rate for bremsstrahlung for the domi- 
nant axial current process is given by ifSTI 

rNN^NNvv ^2nnG^ C\ (3 - cos 0) Ra(oj, q) , (24) 

where n denotes the neutron number density, G the Fermi 
coupling constant, Ca = -ga/2 = -1.26/2 the axial-vector 
coupling for neutrons, and here is the angle between the 
neutrino and antineutrino momenta. The spin or axial re- 
sponse Ra is given by |58 1 
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1 - e-'^/T 



imxA(^, q) , 



(25) 



where Xcr is the spin-density-spin-density response func- 
tion, and we use units with ti - c - - 1 . We use a re- 
laxation time approximation to solve the transport equation 
and obtain the response function [57 1. Under this assump- 
tion, the bremsstrahlung neutrino rates Fi^n^nnvv depend 
on Co-, a quantity connected to the spin relaxation rate Tu- 



Trr 



(26) 



where 



7T m* I 1 



12 ^ 
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^^„^,(k,k')(Tf[(cri +c^2)^ ^^„^,(-k,k')] 



Here ;«* is the effective mass and kp the Fermi momen- 
tum. Also, yicr^.aS)^, k') is the quasiparticle scattering am- 
plitude multiplied by A^(0) = the number of states 
at the Fermi surface, k = pi - p3 and k' = pi - p4 are 
the nucleon momentum transfers, and the average (. . .) is 
over the Fermi surface (for details see ll57l ). The commu- 
tator with the two-body spin operator demonstrates that 
only non-central interactions contribute. The neutrino-pair 
bremsstrahlung and absorption rates /"ata^^a'A'vv are propor- 
tional to Ccr to a good approximation when |w| » I/to-. To 
calculate Co- we treat the strong interaction using the Born 
approximation and evaluate the spin trace in the two-body 
spin space \sms) using a partial-wave expansion. For 5 = 
the spin trace vanishes, thus only 5 = 1 states contribute. 
We derived and expression (see Eq. (6) in Ref. |53|) for 
Co- in terms of partial-wave matrix elements of the strong 
interaction (;?|V^/^|/?), which depend on the relative mo- 
menta p. For all results shown next, we found excellent 
convergence when using 7 partial waves. In addition, we 
checked that the spin relaxation rate for OPE expanded in 
partial waves reproduces the known analytical result (see 
Ref. El). 

In Fig. [16] we show the spin relaxation rate calculated 
from chiral EFT interactions up to N^LO as a function 
of the Fermi momentum over a wide density range {hp - 
1.4fm"' corresponds top = mk^/(3n^) - 1.6xl0''*gcm""^). 
The leading order (LO) contribution includes OPE, which 
provides the standard two-nucleon rates used in current su- 
pernova simulations. The contact terms appearing at LO 
are purely central, thus do not contribute to Co-. Our results 
based on chiral EFT interactions of successively higher or- 
ders from Epelbaum et al. (EGM) |2| show that OPE sig- 
nificantly overestimates Co- for all relevant densities. The 
bands at next-to-leading order (NLO), N^LO, and N^LO 
provide an estimate of the theoretical uncertainty, gener- 
ated by varying the cutoff A as well as a spectral func- 
tion cutoff in the irreducible 27r-exchange ylsp 1 56 2 |. Chi- 
ral EFT interactions at N^LO accurately reproduce low- 
energy NN scattering 12111 . and at this order, Co- is practi- 
cally independent of the N^LO potential for these densities 
(see thin band at N^LO in Fig.[T6]l. Most of the reduction 
of Co- occurs at the NLO level, which includes the leading 
27r-exchange tensor force and shorter range contact terms. 
We observe that at NLO the 2;T-exchange added to the OPE 
enhances the C^, but this is compensated by a further sup- 
pression due to the contact terms. 

We have systematically investigated the sensitivity of 
this observable on different chiral NN forces used in input. 
In the upper panel of Fig. [17] in addition to the band of 
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Fig. 16. Spin relaxation rate given by as a function of Fermi momentum obtained from chiral EFT interactions of successively 
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Fig. 17. The spin relaxation rate given by Co- as a function of Fermi momentum obtained from different chiral N'LO potentials (see 
text) in comparison with the spin relaxation rate obtained from empirical phase shifts (upper panel). The same observable calculated with 
RG-evolved low-momentum interactions obtained from the EM 500 MeV chiral potential is shown in the lower panel for varying 
cutoff A = 1.8 - 2.8 fm" 1. The thickness of the curve for A = 2.0 fm"' is obtained evolving the seven different chiral N^^LO shown in the 
upper panel. The effective mass here is m' = m„ = 939.656 MeV. 



the N-^LO potentials from f2l of Fig. [16] we also show re- 
sults for the Entem and Machleidt (EM) N^LO potentials 
with A = 500 and 600 MeV 1 1 1. The EM 500 MeV poten- 
tial gives results which overlap with the EGM N^'LO band. 



The larger rate obtained with the EM 600 MeV interac- 
tion is due to the failure of the Born approximation in this 
case. We then also used the RG ll60l to evolve N^LO poten- 
tials to low-momentum interactions Viowk with A = 1.8- 
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2.8 fm % to extend the cutoff variation and estimate of the 
theoretical uncertainty. The results are shown in the lower 
panel of Fig. [T7] The RG evolution preserves the long- 
range pion exchanges and includes sub-leading contact in- 
teractions, so that NN scattering data are reproduced. The 
thickness of the curve for A = 2.0 fm"' is obtained af- 
ter evolving all seven N-^LO potentials of the upper panel. 
This shows the universality of Viow* and that, after high- 
momentum modes are integrated out, the particle-particle 
channel (for EM 600 MeV) is rendered more perturba- 
tive [60 1 . 

At low densities, two-nucleon collisions dominate and 
NN scattering observables provide a model-independent 
result for the spin relaxation rate. To study this we can re- 
place the strong potential V by the T matrix, which is pa- 
rameterized in terms of empirical phase shifts and mixing 
angles lIMll . taken from the Nijmegen Partial Wave Analy- 
sis (nn-online.org). For all densities in Fig. [17] we find that 
the rates obtained in the Born approximation from chiral 
N^LO potentials and V\ovik interactions are close to those 
from phase shifts. Combined with the RG insights, this 
demonstrates that the non-central part of the strong neu- 
tron - neutron amplitude is perturbative for these lower 
cutoff interactions. 

The relatively small sensitivity of Co- to the NN poten- 
tial used in input, allows us to fit a very simple analytical 
function 

Ca_ ^ 0.86(^F/fm-')3-^ 
MeV 1 -H 10.9(;tF/fm"')3 f> ' 

that can be used in hydrodynamic simulations of super- 
novae explosions. In fact, our main finding is that two- 
pion exchange interactions and shorter-range non-central 
forces, included beyond the LO in EFT, reduce the neu- 
trino rates significantly, with respect to the standard rates 
used in supernovae simulations. Future work will include 
neutron-proton mixtures and charged currents, to system- 
atically improve the neutrino physics input for astrophysics. 
We also plan to investigate other many-body eff'ects on 
collisions including 3NFs, which are expected to be im- 
portant at higher densities. 



5 Conclusions and Outlook 

The investigation of electroweak reactions involving light 
nuclei has entered a new era, due to new theoretical ad- 
vances and to the number of physical aspects one can study. 
It allows one to tackle a variety of questions, ranging from 
more fundamental ones on how nucleons interact with each 
other and spanning to astrophysics, where one can learn 
about the role of nuclear input in astrophysical objects. 

Here, we have reported on the recent progress done 
in calculations of electroweak reactions on light nuclei, 
mainly using the LIT method in conjunction with EIHH 
expansions of the wave function. In such way, bound and 
scattering observables are treated on an equal footing. The 
use of the LIT method is so far the only viable procedure 
to access information on the inelastic response of A > 3 
nuclei for energies beyond the three-body disintegration 



threshold. A direct calculation of all involved channels is 
in fact impossible. The LIT method, by circumventing the 
problem via an integral transform and turning it to the so- 
lution of a bound state equation, constitutes a very practi- 
cal and yet effective approach to obtain the exact response 
function for light nuclei. 

The excitation of nuclei with external electromagnetic 
probes is a fundamental tool to get information about the 
nuclear dynamics via a comparison of theoretical predic- 
tions with experiments. Theoretical efforts need to be ad- 
dressed to the identification of many-body observables that 
are sensitive to 3NFs, that are not constrained by NN scat- 
tering. We have shown that the photo-absorption reaction 
on '*He presents a relatively mild sensitivity to 3NF in the 
energy range of the cross section peak. Furthermore, the 
experimental situation on the a-particle does not allow one 
to discriminate among difference force models. On the other 
hand, the electron scattering reactions can potentially pro- 
vide us with useful information on the nuclear dynamics. 
Here in fact, we find a larger effect of 3NF, especially 
at low-momentum transfer. In this region no experimen- 
tal measurements are found in literature, but new data have 
been taken in Mainz and are to be analyzed ll35l . More 
work needs to be done from the theoretical point of view 
by using other potential models available and by perform- 
ing further sensitivity studies. For example, calculations of 
Ri and Rj at low q with chiral EFT potentials are called 
for. If the new experimental data are precise enough, one 
could explore the possibility of calibrating low energy con- 
stants on such continuum observables. They are in fact ex- 
pected to contain a richer information about nuclear dy- 
namics than simple bound stable observables, like the bind- 
ing energy, due to the fact that many channels are involved. 

We have also shown recent results for weak observ- 
ables where excitations of nucleonic matter are induced 
by neutrino interactions. Such processes are very impor- 
tant in astrophysical scenarios, like in supernovae explo- 
sions. We have reported on a recent ab-initio calculation 
li46J of the v-a scattering with two and three-body forces, 
that can be used as input for supernovae simulations. One 
can observe that such astrophysically relevant observables 
also display a dependence on three-body forces, which is 
interestingly similar to that found in electron scattering. 
We also have shown how a recent study [.49 J pointed out 
the importance of A = 3 disintegration as a way to con- 
tribute to the neutrino energy transfer to the shock wave 
in supernovae. Finally, we have presented a study of the 
neutrino-bremsstrahlung process using chiral EFT poten- 
tials. We have shown how modern two-body Hamiltonians 
can improve our understanding of such reaction rates. Our 
results can be incorporated in simulations and are part of a 
larger project to systematically improve the nuclear micro- 
scopic input for supernovae simulations. Future work will 
address the role of many-body physics, including the effect 
of 3NFs, relevant at higher densities. Summarizing, since 
often no experiments are possible to simulate astrophysical 
conditions, all presented calculations show how few-body 
physics can provide valuable input for hydrodynamic sim- 
ulations of supernovae explosions. 
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